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\Q [ Abstract 

Using the Bayesian theory of model comparison, a new cosmological 
CN ! model due to John and Joseph [M. V. John and K. Babu Joseph, Phys. 

Rev. D 61 87304 (2000)] is compared with the standard Q\ ^ cosmolog- 
t— i ■ ical model. Their analysis based on the recent apparent magnitude-redshift 

data of Type la supernovas found evidence against the new model; our more 
Q careful analysis finds instead that this evidence is not strong. On the other 

\ hand, we find that the angular size-redshift data from compact (milliarcsec- 

Q-V ond) radio sources do not discriminate between the two models. Our analysis 

O ' serves as an example of how to compare the relative merits of cosmological 

+3 . models in general, using the Bayesian approach. 

PACS No(s): 98.80.Es, 02.50.-r 
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I. INTRODUCTION 

In a recent publication QTJ , it was argued, by modifying an earlier ansatz 
by Chen and Wu 0, that the total energy density p for the universe should 
vary as a~ 2 where a is the scale factor of its expansion. If the total pressure 
is p, then this argument leads to p + 3p = for the universe. This deduction 
was made possible by the use of some dimensional considerations in line 
with quantum cosmology. The reasoning is as follows: Taking the comoving 
coordinate grid as dimensionless, we attribute a distance dimension to the 
scale factor a. Since there is no other fundamental energy scale available, one 
can always write p as Planck density (p p i = c 5 /hG 2 = 5.158 x 10 93 g cm~ 3 ) 
times a dimensionless product of quantities. The variation of p with a can 
now be written as 

hi 

p oc p p i — , 
[a 

where l pl = (frG/c 3 ) 1 / 2 = 1.616 x 1(T 33 cm is the Planck length. It is easy to 
see that n < 2 (n > 2) will lead to a negative (positive) power of h appearing 
explicitly on the right hand side of the above equation. It was pointed out 
that such an ^-dependent total energy density would be quite unnatural in 
the classical Einstein equation for cosmology, much later than the Planck 
time. However, the case n = 2 is just right to survive the semi classical limit 
H — > 0. Thus it was argued that if we take quantum cosmology seriously, 
then p oc a~ 2 or equivalently p + 3p = 0, for a conserved p. Solving the 
Friedmann equations gives a coasting evolution for the universe; i.e., 

a = mt, 

where m = \Jk/(Vt — 1) is a proportionality constant; Vt is the total density 
parameter and k — 0, ±1 is the spatial curvature constant. 

It shall be noted that p + Sp = is an equation of state appropriate for 
strings or textures and that it is unrealistic to consider the present universe as 
string-dominated. But in [[[]], it was shown that this ansatz will lead to a real- 
istic cosmology if we consider that p comprises of more than one component, 
say, ordinary matter (relativistic or nonrelativistic) with equation of state 
p m = w p m and a cosmological constant A, which is time-varying. Let p\ 
denote the energy density arising from A and p\ = —p\ be the corresponding 
pressure. With 



3 



P = Pm + PA, P = Pm+PA, 

the condition p + 3p = will give 

Pm _ 2 

Pa 1 + 3w ' 

and this gives a realistic model for the universe. It was also shown that 
this simplest cosmological model is devoid of the problems like the horizon, 
flatness, monopole, cosmological constant, size, age of the universe and the 
generation of density perturbations on scales well above the present Hubble 
radius in the pure classical epoch. The solution of the cosmological constant, 
age and density perturbation problems deserve special mention since these 
are not solvable in an inflationary scenario. Moreover, the evolution of tem- 
perature in the model is nearly the same as that in the standard big bang 
model and if we assume the values fl m = 4/3 and fl\ = 2/3, then there is 
no variation in the freezing temperature with the latter model, and this will 
enable nucleosynthesis to proceed in an almost identical manner. It also may 
be noted that an almost similar model which predicts the above values for 
the density parameters was proposed earlier 0, from some more fundamental 
assumptions based on entirely different grounds. 

However, it should be remarked that the argument given above, which 
leads to this cosmology, is heuristic and not based on formal reasoning. It 
should be taken only as a guiding principle. Also we note that it has some 
unusual consequences like the necessity of continuous creation of matter from 
vacuum energy, though it was argued in [|1], that such creation will be too 
inaccessible to observation. 

But it was mentioned in [|l| that, in spite of the those successes in predict- 
ing observed values, the recent observations of the magnitudes of 42 high- 
redshift Type la supernovas |§] is a set back for the model. A statement was 
explicitly made to the effect that the predictions of Q m and fl\ for the present 
model are outside the error ellipses given in the fl m — plot in |{|] and it 
was claimed that this discrepancy is a serious problem. In this paper, we 
study this issue in detail to see how strong is the evidence against this model 
when compared with the standard model with a constant A ^ 0, discussed 
in f|, ||. Jackson and Dodgson || |7j have examined the latter model in the 
light of Kellerman's || and Gurvits' ||] compilations of angular size-redshift 
data for ultracompact (milliarcsecond) radio sources. Gurvits' compilation of 
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such data, which are measured by very long-baseline interferometry (VLBI), 
is claimed to have no evolution with cosmic epoch. Several authors (for e.g., 
flCf ) have made use of these data to test their cosmological models. In the 
present paper, we also analyze the Gurvits' data to test the new model. Us- 
ing the Bayesian theory of statistics, we compare the new model discussed 
above with the standard model with a non-zero cosmological constant, us- 
ing both the apparent magnitude-redshift data and the angular size-redshift 
data. It is found that there is no strong evidence against the new model when 
the apparent magnitude-redshift data are considered. This is contradictory 
to the statement made in The angular size-redshift data, on the other 
hand, are found to provide equal preference to the standard model and the 
new one. 

The remainder of this paper takes the new theory as given and compares it 
with other standard cosmological models. The analysis shall be viewed as an 
example of using Bayesian theory to test the relative merits of cosmological 
models, a method which is claimed to have many positive features when 
compared to indirect arguments using parameter estimates. As such, the 
technique described here has wider applicability than just to the comparison 
of two cosmological models. 

The paper is organized as follows. In Section II, we discuss the Bayesian 
theory of model comparison for the general case. Section III discusses com- 
parison of the two models using apparent magnitude-redshift data and in 
Section IV, we compare the models with the angular size-redshift data. Sec- 
tion V comprises discussion of the results. 

II. BAYESIAN THEORY OF MODEL COMPARISON 



The Bayesian theory of statistics [11, 12] is historically the original ap 



proach to statistics, developed by great mathematicians like Gauss, Bayes, 
Laplace, Bernoulli etc., and has several advantages over the currently used 
long-run relative frequency (frequentist) approach to statistics, especially in 
problems like those in astrophysics, where the notion of a statistical ensem- 
ble is highly contrived. The frequentist definition of probability can only 
describe the probability of a true random variable, which can take on vari- 
ous values throughout an ensemble or a series of repeated experiments. In 
astrophysical and similar problems, ensembles and repeated experiments are 
rarely possible and we speak about the probability of a hypothesis, which 
can only be either true or false, and hence is not a random variable. The 
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Bayesian theory will help assign probabilities for such hypotheses by consid- 
ering the (often incomplete) data available with us. For example, Laplace 
used Bayesian theory to estimate the masses of planets from astronomical 
data, and to quantify the uncertainty of the masses due to observational er- 
rors [l3|]. In fact, this theory finds application in all those problems where 
one can only have a numerical encoding of one's state of knowledge. 

In the Bayesian theory of model comparison, it is common to report model 
probabilities via odds, the ratios of probabilities of the models. The posterior 
(i.e., after consideration of the data) odds for the model Mj over Mj are 

Uij P (M 3 \D,iy 

where p(Mi\D,I) refers to the posterior probability for the model Mj, given 
the data D and assuming that any other information I regarding the models 
under consideration is true. Using Bayes's theorem, one can write the above 
equation as 

p(Mi\I)C(Mj) 
Uij p{M j \I)C{M j y [L) 

where p(Mi\I) is called the prior probability; i.e., any probability assigned to 
the model Mj before consideration of the data, but assuming the information 
I to be true. When / does not give any preference to one model over the 
other, these prior probabilities are equal so that 

° ii = l^)= Bii - (2) 

By is called the Bayes factor. £(Mj) denotes the probability p(D\Mi) to 
obtain the data D if the model Mj is the true one and is referred to as the 
likelihood for the model Mj. The models under consideration will usually 
have one or more free parameters (like the density parameters f2 m , etc. 
in the case of cosmological models), which we denote as a, (3, .. . £(Mj) can 
be evaluated for models with one parameter as 

£(M) = p(D\Mi) = J da p(a|M)A(a), (3) 

where p(a\Mi) is the prior probability for the parameter a, assuming the 
model Mj to be true. Ci(a) is the likelihood for a in the model and is 
usually taken to have the form 
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A(a) = exp[-xV)/2]. (4) 

where 

k \ a k J 

is the x 2 statistic. Here A k are the measured values of the observable A, 
A k (a) are its expected values (from theory) and a k are the uncertainties in 
the measurements of the observable. 

Generalization to the case of more than one parameter is straight forward. 
As a specific case, consider a model Mi with two parameters, a and (3, having 
flat prior probabilities; i.e., we assume to have no prior information regarding 
a and f3 except that they lie in some range [a, a + Aa] and \f3, (3 + A/3], 
respectively. Then p(a\Mi) = 1/Aa, p((3\Mi) = l/A/3 and hence 

X X f a+Aa rP+A/3 

£(Mi) = A^ L d ^M-X 2 {a^)/2\. (6) 

It is instructive to rewrite this equation as 



C(Mi) = — / dad 



a 



In this case, 



1 rP+W 

Cl{a) = A3J d(3exp[- X 2 (a,(3)/2} 



A/3 J p 

is called the marginal likelihood for the parameter a. 

A. Interpretation of the Bayes factor 

The interpretation of the Bayes factor B^^ which is given by @ and 
which evaluates the relative merits of model Mj over model Mj, is as follows 



14| : If 1 < Bij < 3, there is an evidence against Mj when compared with Mj, 



but it is not worth more than a bare mention. If 3 < B^ < 20, the evidence 
against Mj is definite but not strong. For 20 < B^ < 150, this evidence is 
strong and for By > 150, it is very strong. 



III. COMPARISON USING REDSHIFT-MAGNITUDE DATA 
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For an FRW model which contains matter and a cosmological constant, 
the likelihood for these parameters, i.e., £i(Q m ,QA) can be assigned using 
the redshift-apparent magnitude data in the following manner |14]]. Before 
consideration of the data, let us agree that fl m lies somewhere in the range 
< Q m < 3, Qa in the range — 3 < Q\ < 3 and accept this as the only 
prior information /. Let fi k be the observed best-fit distance modulus for the 
supernova number k, s k its uncertainty and z k is the cosmological redshift, 
with w k its uncertainty. We can write the expression for \ 2 as 



X 




— ^k 
O'k 



Here, 



with 



h = Hk + n k = g k - rj + n k , 



(7) 



(8) 



l^k = 9k-V = ^ log 



Dl(z; ^A, Hq) 



+ 25 



IMpc 

being the redshift-apparent magnitude relation. The luminosity distance is 

where c is the velocity of light, H 



DAz\ Urn, $l\,Hn) = cHn l d Tl {z\ fi m , fi/ 



is the Hubble constant at the present epoch and dt is the dimensionless 
luminosity distance. g k = g(z k ) is the part of [i k which depends implicitly 
on VL m and Vt^ and 77 is its ifo-dependent part. The latter quantity can be 
written as rj = 51og(/i/c2) — 25 where H Q = h x 100 km s -1 Mpc" 1 and C2 
is the speed of light in units of 100 km s -1 . The probability distribution 
for the value n k in Equation (||) is assumed to be a zero-mean Gaussian 
with standard deviation a k , where a k = si + \z k )] 2 w kl in the absence of 
systematic or evolutionary effects. 

One can evaluate C(Q m , Q\, rj) in a manner similar to that in @, where 
X 2 now is a function of the three parameters fl m , fl\ and rj. The likelihood 
for fl m and fl\, denoted as £(Q m ,QA) can be obtained by the technique of 
marginalising over r], if one assumes a flat prior probability for r] in some 
appropriate range. 

To do this, we define = sjY^kO-l a k) where s is the posterior uncer- 
tainty for rj and let I/A7/ a flat prior probability be assigned to rj. (These 
being the same for all models, will get canceled when evaluating probability 
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ratios.) Using these definitions, the marginal likelihood (defined at the end 
of Sec. II) for the density parameters is 

£(n m ,n A ) = ^-Jd V e-* 2 / 2 . (9) 

Evaluating this integral analytically flyj], one assigns a likelihood for the 
parameters Q m and fl\ in any one model as 

C(Q m ,n A ) = ^^e-^ 2 , (10) 



where 



g(n m ,n A ) = 2^ 2 ' ( n ) 

is of the form of a x 2 -statistic, with f) the best fit (most probable) value of 
77, given fl m and The latter can be computed as jl4| 

k a k 

Now, we compare the model in [|, ||] (model M 1? having parameters fl m , 
Q A and rf) with the new model discussed in Sec. I (model M 2 , having only 
the parameters Q m and rf). The Bayes factor B 12 can be written with the 
help of Eq. (|) and Eq. (|) as 

= C(M 1 ) _ JdQ m J^ap^^aIMQA^^a) 
12 £(M 2 ) J^ m p(l] m |M 2 )£ 2 (O w ) 1 j 

With the information I at hand, one can assign flat prior probabilities 
j9(fi m ,fiA|M!) = 1/18 and p{Vt m \M 2 ) = 1/3 . Using Eqs. (|) and (|TTJ) we 
can write the above as 

_ /-3 rf ^A Jo exp[-gi(fl m , Oa)/2] 
6/ 3 rffi m exp[-g 2 (fi m )/2] 

Our first step in the evaluation of B i2 is to find q given in Eq. (|IT|), for 
both the models. For Model 1, we have to use 

(/(*) = 5 log{(l + z) l^r^sinnO^I 1 / 2 /^)]}, 
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where fl k = 1— and sinn(a;) = sinx for Q m +Q A > 1, sinn(x) = sinhx 
for Q m + Q\ < 1 and sinn(x) = x for Q m + Q\ = 1. Also 

I(z) = f [(1 + z') 2 (l + ®< m z') - z'{2 + z')(^a)]~ 1/2 dz'. 
Jo 

For Model 2, the function g(z) can be written as 

1 

(7(2) = 5 log{m(l + z)sinn[ — ln(l + z)}}, 

where m = \jlkj (3Q m — 2) for the nonrelativistic era and sinn(x) = sinx 
for Q m > 2/3, sinn(x) = sinhx for Q m < 2/3 and sinn(a;) = x for fl m = 2/3. 

Using these expressions, Eq. (0) is numerically evaluated to obtain B u = 
3.1. (In this calculation, we have used the data corresponding to the Fit C 
in 0, which involve 54 supernovas.) As per the interpretation of given 
in Sec. II. A, the above is an evidence against Model 2, but it is only barely 
definite; the discrepancy is not a "serious problem" as had been stated in |T|. 

IV. COMPARISON USING ANGULAR SIZE-REDSHIFT DATA 

For this purpose, we use the Gurvits' data and divide the sample which 
contains 256 sources into 16 redshift bins, as done by Jackson and Dodgson 
and shown in their Fig. 1 0. For Model 1, we use the expression for angular 
size 



&9 = — = i 

d A (1 + z )-i (k/n k ) 1/2 ^ S inn[|^|V2/'(^)] 
dH (1 + z) 



where 



c (k/n k ) 1/2 smn[\n k \V>F(z)} 



(15) 



"l+z d T 
n*) = l —TTJ2- (16) 



x 



(p k + Vt m x + 



x 2 J 



Here d is the linear dimension of an object, d A is the angular size distance 
and and sinn(a;) are defined as in the case of Model 1 in the last section. 
Similarly for Model 2, we have 
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M= d_ = dH (1 + z) 

d A cm sinnfi ln(l + z)\ ' 

where m and sinn(x) are defined as in the earlier case of Model 2. For the 
purpose of comparison, we only need to combine the unknown parameters 
d and H to form a single parameter p = dH /c. Thus Model 1 has three 
parameters p, Q m and whereas Model 2 has only the parameters p and 
Q m . As in the previous case, we accept < Q m < 3 and — 3 < < 3 as the 
prior information /. With these ranges of values of Q m and P is found to 
give significantly low values of \ 2 on ly f° r the range 0.1 < p < 1 in both the 
models, p being given in units of milliarcseconds. The formal expressions to 
be used are 




M r/i (is) 

and 



_ £(Mi) ^ aL An A / dpJdVl m J dn A exp[-xj(p, gm, ^a)] 
12 £(M 2 ) ^^/dp/dfi m exp[-xl(p,fi m )] 

6 Jo 1 .! rf P /o 3 exp[-x|/2] 

The result obtained is Bu ~ 1. This may be interpreted as providing 
equal preference to both models. 



V. DISCUSSION 



While evaluating the Bayes factors using both kinds of data, we have 
assumed that our prior information / regarding the density parameters is 
< Q m < 3 and — 3 < ^a < 3. The range of values of Q\ considered in 
@| is -1.5 < fi A < 3 and in it is -4 < tt A < 1. Even if we modify the 
range of this parameter in our analysis to some reasonable extent, the main 
conclusions of the paper will remain unaltered. For example, if we accept 
< Q m < 3 and —1.5 < ^a < 1-5 as some prior information I', the Bayes 
factors in each case become 3.8 and 0.8, in place of 3.1 and 1, respectively. 
Instead, if we choose I' as < Q m < 3 and — 6 < ^a < 6, the corresponding 
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values are 1.55 and 1.4, respectively. These do not change our conclusions 
very much in the light of the discriminatory inequalities mentioned in Sec. 
II.A. 

In order to get an intuitive feeling why the standard (Mi) and new 
(M 2 ) models have comparable likelihoods, consider Figs. [I] and ||. Fig. 
|l] is for the apparent magnitude-redshift data and plots the quantities L' = 
| (if^A exp [— qi{yt m) Vt^)/2] (curve labeled Mi) and L' = exp[— q 2 {Vl m )/2] 
(curve labeled M 2 ) against Q m . From the definition of marginal likelihood 
given at the end of Sec. II and from Eqs. (9)-(14), it can be seen that these 
two curves correspond to the marginal likelihoods for the parameter Vt m in 
models Mi and M 2 , respectively (apart from some multiplicative constants, 
which cancel on taking ratios). Similarly, Fig. [2|, which is for the angular 
size-redshift data, plots C = 6x q 9 Jo 1 ,! dp f_ 3 d£l\ exp[— (curve Mi) and 
C = Jo ! dpexp[— x|/2] (curve M 2 ) against Vt m . Eq. (|19|) allows us to in- 
terpret these terms as the marginal likelihoods for Q m in models Mi and M 2 , 
respectively. In fact, these curves rigorously show the integrands one must 
integrate over fl m to get the Bayes factors. Using the apparent magnitude- 
redshift data, a lower value of value of q (which is a modified x 2 statistic) 
is obtained for model Mi whereas for angular size-red shift data, lower \ 2 is 
claimed by model M 2 . However, the areas under the curves are comparable 
in both cases and this shows why the Bayes factors are also comparable. This 
is one of the strong points of the Bayesian method, in contrast to frequentist 
goodness of fit tests, which consider only the best fit parameter values for 



comparing models |Tl 



These figures, however, show some feature that is disturbing for the new 
model. Figs. [1] and || indicate best fit values of Q m = and fl m = 0.42, 
respectively, for this model. In both cases it appears to rule out the value 
Q m = 4/3 that is needed to meet the constraints on nucleosynthesis, a con- 
dition which had been stated in the introduction. Though, as mentioned 
above, Bayesian model comparison does not hinge upon the best fit values in 
evaluating relative merits of models, one would desire to have an agreement 
between predicted and observed parameter values. A natural option in such 
cases would be to compare the models by adjusting the prior regarding the 
parameters so that any additional information is accounted for. But we have 
not attempted this in our analysis. 

The constant Oa ^ model we considered has one parameter in excess 
of the new model in both cases. It should be kept in mind that in the 
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Bayesian method, simpler models with less number of parameters are often 
favored unless the data are truly difficult to account for with such models. 
Bayes's factors thus implement a kind of automatic and objective Occam's 
razor. In this context, it is interesting to check how the new model fares 
when compared with flat (inflationary) models where O m + ^A = 1, by which 
condition the number of parameters of model Mi are reduced by one. This 
makes the two models at par with each other, with regard to the number 
of parameters. We have calculated the Bayes factor between this flat model 
Mi and the new model M 2 , using the apparent magnitude-redshift data and 
the result is Bi 2 = 5.0. This appears to be a slightly more definite evidence 
against the new model than the corresponding result obtained in Sec. Ill 
(B\2 = 3.1). (However, inflationary models with a constant A-term suffer 
from the 'graceful exit problem' for A; i.e., in order to explain how A manages 
to change from its GUT magnitude to ~ 10~ 126 of its initial value, some 
extreme fine tuning would be required |15|). On the other hand, a comparison 
of Q m + Q\ = 1 model with the new model using angular size-redshift data 
gives a value for the Bayes factor B 2 i = 15, which shows that this data is 
difficult to account with the flat inflationary models than with the new one. 
The results we obtained, while using the information I, are summarized in 
Table I. 

When compared to the frequentist goodness of fit test of models, which 
judges the relative merits of the models using the lowest value of x 2 (even 
when it is obtained by some fine tuning or by having more parameters), the 
present approach has the advantage that it evaluates the overall performance 
of the models under consideration. The Bayesian method is thus a very 
powerful tool of model comparison and it is high time that the method is 
used to evaluate the plausibility of cosmological models cropping up in the 
literature. It is true that since we have only one universe, one can only resort 
to model making and then to comparing their predictions with observations. 
Again, since we cannot experiment with the universe, it is not meaningful to 
use the frequentist approach. We believe that the only meaningful way is to 
use the Bayesian approach in such cases. Here we have made a comparison 
between the model in 0, || with the new model in |l|. It deserves to be 
stressed that the recent apparent magnitude-redshift observations on Type 
la supernovas do not pose a "serious problem" to the new model, as had been 
claimed in 0. The angular size-redshift data, on the other hand, do not 
discriminate between the general Q\ ^ model and the new model and they 
provide definite but not strong evidence against standard flat (Q m + Qh = 1) 
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model when compared to the new one. 

Here it is essential to point out that Bayesian inference summaries the 
weight of evidence by the full posterior odds and not just by the Bayes fac- 
tor. Throughout our analysis above, we have assumed that the only prior 
information with us is either I (stated in the beginning of Sec. Ill) or V 
(stated in the beginning of Sec. V), which helps to make the posterior odds 
equal to the Bayes factor. However, when the Bayes factor is near unity, the 
prior odds p(Mj | I)/p(Mj \ I) in Eq. (|l|) become very important. The stan- 
dard f^A 7^ model and the standard flat (inflationary) models are plagued 
by the large number of cosmological problems (as mentioned in Sec. I) and 
the new model has the heuristic nature of its derivation and the problem 
with nucleosynthesis, setting (subjective) prior odds against each of them. 
In the context of having obtained comparable values for the Bayes factor, the 
Bayesian model comparison forces us to conclude, in a similar tone as in ||14j| , 
that the existing apparent magnitude or angular size-redshift data alone are 
not very discriminating about these cosmological models. It is also worth re- 
marking here that the Bayesian theory tells us how to adjust our plausibility 
assessments when our state of knowledge regarding an hypothesis changes 



through the acquisition of new data 11 1. Concerning future observations, 



one would have to say that if the supernova test is extended to higher red- 
shifts and if the astronomers are sure about the standard candle hypothesis, 
then the theories can be tested for such new data using Bayesian model com- 
parison, using what we have now obtained as the prior odds. In this context, 
it also deserves serious consideration to extend the analysis made here to 
other cosmological data, like that of cosmic microwave background radia- 
tion and primordial nucleosynthesis. Hopefully, further analysis and future 
observations may help to give more decisive answers on these questions. 
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Table I. 



Data Model Mi Model M 2 Bayes factor Interpretation 



m — z 


Standard 

Oa ^0 model 


New model 


B\2 


= 3.1 


Slightly definite but 
not strong evidence 
against the new model 


m — z 


Standard flat 
Oa ^0 model 


New model 


B\2 


= 5 


Definite but 

not strong evidence 

against the new model 


6-z 


Standard 

Oa t^O model 


New model 




= 1 
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Figure 1: V vs O m for both models, using the apparent magnitude-redshift 
data for Type la supernova. The curves Mi and M 2 correspond to the 
marginal likelihoods for fl m for the standard fl\ ^ model and the new 
model, respectively (apart from some multiplicative constants, which cancel 
on taking ratios). 



REFERENCES 



17 




Figure 2: Marginal likelihood vs VL m for both models, using the angular size- 
redshift data. The curves M\ and M 2 correspond to the marginal likelihoods 
for Q m for the standard fl\ ^ model and the new model, respectively. 



